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Abstract 

We present a detailed study of the cohider observable mT2 applied for pair-produced 
superparticles decaying to visible particles and a pair of invisible lightest supersym- 
metric particles (LSPs). Analytic expressions of the maximum of mT2 over all events 
(m™|^) are derived. It is noticed that if the decay product of each superparticle 
involves more than one visible particle, m™!^ being a function of the trial LSP mass 
has a kink structure at = true LSP mass, which can be used to determine 
the mother superparticle mass and the LSP mass simultaneously. To see how well 
m^2^ can be constructed from collider data, a Monte-Carlo analysis of the gluino 
mT2 is performed for some superparticle spectra. 



1 Introduction 



The Large Hadron Collider (LHC) at CERN will explore soon the TeV energy scale where 
new physics beyond the Standard Model (SM) is likely to reveal itself [HIS]. Among various 
proposals, weak scale supersymmetry (SUSY) [3] is perhaps the most promising candidate 
for new physics at TeV as it provides a solution to the gauge hierarchy problem while 
complying with gauge coupling unification. Furthermore, with R-parity conservation, the 
lightest supersymmetric particle (LSP) becomes a natural candidate for the non-baryonic 
dark matter in the Universe. 

Once SUSY signals are discovered through event excess beyond the SM backgrounds 
in inclusive search channels, the next step will be the measurement of SUSY particle 
masses and other physical properties through various exclusive decay chains [1] . Then, it 
might be possible to reconstruct the underlying SUSY theory, in particular the soft SUSY 
breaking terms, using the observed SUSY particle masses. In this regard, experimental 
information on gaugino masses can be particularly useful for distinguishing different SUSY 
breaking schemes as the theoretical predictions of low energy gaugino masses are quite 
robust compared to those on sfermion masses [5]. 

In a recent paper [6], we have examined the collider observable "gluino 771^2 (strans- 
verse mass)" which corresponds to the Cambridge mT2 variable [131 [H] applied to pair 
produced gluinos each of which is decaying to two quarks and one invisible neutralino 
LSP. Analytic expression of m™2^ (= maximum of the gluino mT2 over all events) as a 
function of the trial LSP mass has been discussed with an observation that m™2^ has 
a kink structure, i.e. a continuous but not differentiable cusp, at = true LSP mass, 
from which the gluino mass and the LSP mass can be determined simultaneously with 
good accuracy. If squarks are lighter than gluino, the gluino m™2^ could determine the 
squark mass also. In this paper, we wish to provide a detailed discussion of mT2 in more 
general context, including the features which have been reported in [6j*l . 

This paper is organized as follows. In section 2, we discuss some generic properties 
of mT2, and derive the analytic expression of m™|'^ for general symmetric decay of pair- 
produced mother superparticles. In section 3, we consider two specific processes, the 
decays of squark pair and of gluino pair, to examine the structure of mT2 in somewhat 
detail, and also perform a Monte Carlo LHC simulation for some superparticle spectra to 
examine how well can mT2 and m™|^ be constructed from real collider data. Section 4 is 
devoted to the conclusion. 

* After [6], the kink structure of the endpoint values of transverse mass observable has been discussed 
also in [7l[8]. For other approaches to measure superparticle masses at hadron collider, see [4} l9l [TOl [TT|. 
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2 Generic features of mT2 



The kinematic variable 'transverse mass' (mj^) has been introduced to measure the W 
boson mass from the decay W ^ lu [12], for which the transverse mass is given by 

ml = m1 + ml + 2{El^E^ - p?^ ■ p^), (1) 

where m^, rriy and p^, p^ denote the mass and transverse momentum of the corresponding 
particle, respectively, and the transverse energies are defined as 

Elr=^J\v'T? + ml El^ = ^\^'^^\^ + ml. (2) 
On the other hand, the physical W mass is given by 

= mf + ml + 2pi ■ 

= rri^ + ml + 2(E^E^ cosh Ar/ - p^ ■ p^) > m^, (3) 

where A?7 = rji — rji, is the rapidity difference for a 4-momentum parameterized as 

= {Et cosh ?7, Pt, Et sinh 77) . 

Although the neutrino cannot be observed directly, its transverse momentum can be 
inferred from the measured total missing transverse momentum in event-by-event basis. 
Then, for each event, the transverse mass of W can be constructed from the observed 
values of p^ and p^, and an endpoint measurement of the ijit distribution determines 
the mother particle mass mw'- 

mw = max [ ] . (4) 
{all events} 

More challenging situation for experimental measurement of unknown mother particle 
mass would be the case when there are more than one final state particles escaping the 
detection, so that the transverse momentum of each invisible particle can not be deter- 
mined although the total missing transverse momentum is known. Additional difficulty 
would arise if the mass of the invisible daughter particle were not known in advance. 
Such situation is what one actually encounters in supersymmetric extension of the SM 
with conserved -R-parity, in which superparticles are pair produced in collider experiment 
and each superparticle decay ends up with producing an invisible lightest supersymmetric 
particle (LSP) with unknown mass. 

The 771x2 variable [131 [H] which is sometimes called the "stransverse mass" is a gen- 
eralization of the transverse mass to the case that a pair of massive mother particles are 
produced in hadron collider with a vanishing total transverse momentum in the laboratory 
frame, and subsequently decay to daughter particles including two invisible particles in 
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the final statqll. In this paper, we will concentrate on the case that each mother particle 
decays into the same set of daughter particles, since such symmetric decay typically has 
higher event rate while showing the non-trivial structure which will be discussed in the 
following. Fig. 1 shows an example of such process in which mother superparticles were 
pair-produced and each of them decays into one neutralino LSP (x?) and some visible 
particles. While the invisible part of each decay consists of only one particle (neutralino 
LSP), the visible part might contain one or more visible particle(s) in general. 




Figure 1: Kinematic situation for mT2 where p^*^* denotes the total missing transverse 
momentum. 

With two invisible LSPs in the final state, each LSP momentum can not be determined 
although the total missing transverse momentum p™^* can be measured experimentally. 
Furthermore, the LSP mass might not be known in advance. In such situation, one can 
introduce a trial LSP mass m^, and define the mT2 variable as follows [13], [1] 

/ visil) (1) vis(2) (2) X 

mT2{PT , Pr % ^viU ^x) = mm 



maxjm^'*, m^'*} , (5) 



where the minimization is performed over trial LSP momenta p^^*'' constrained as 

„X(1) I „X(2) _ miss 

Pt "t Pt — Pt i 



^In Ref.|14j. mT2 has been further generalized to the case involving more missing particles than two. 



3 



and m^'' [i = 1,2) denotes the transverse mass of the decay product of each initial mother 
particle: 



m« = ^J{m^:ly + + 2( E^r^'^^ Ef^ - p^^^^ ■ P^^'^), (6) 

where rrS^J^ and p^*'*'-*^ are the total invariant mass and the total transverse momentum 
of the visible part of each decay product: 

vis{i) _ \ ^ 
Pt — / ,PaT, 

a 

{^^1? = $^m^ + 2^(i?„i?^-p„-p^), (7) 

a a>f3 

where E^ = \/ + Ip^p, Pa, and PaT denote the energy, momentum, and transverse 
momentum, respectively, of the a-th visible particle in the decay product of one mother 
particle, which are measured in the laboratory frame. Here, we ignore the initial state 
radiation effect, and then the total missing transverse momentum is given by 

p-- = -(p«^W+p«^(2)^^ 

and the transverse energies of the each visible system and of the LSP are defined as 



In fact, one can consider also an 712x2 with m^*^^ replaced by the transverse mass of the 
visible part: 

(m^^^W)^ = Y,< + 2Y,{EaTEpT - VaT ■ P/3t), (10) 

a a>f} 

where E^t = \/ m"^ + Ip^tP- Since py*^*-* and m|,*g are treated as independent variables 
in the definition ([5]), once one finds the functional form of 771^2 in terms of p^*'**'*'' and m^*^, 
the other mT2 defined in terms of p™*^*'* and m™*''*^ can be easily obtained by replacing 

(*) -4-1, vis(i) 

m^-g with rrij, 

Even without knowing the transverse momentum of each LSP, one can ensure that 
the true mother particle mass rh can not be smaller than when the trial LSP mass 
is chosen to be the true LSP mass m^o. This suggests that rh might be able to be 
determined by the endpoint value of 171x2 distribution: 



max/ \ / vis(l) (1) vis(2) (2) \ 

{ail events} 
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which is a function of the trial LSP mass m^, satisfying 

m™2^(m^ = m^o) = m = mother particle mass. (12) 

In order to see how 771^2 is determined for a given event, let us first consider the 
minimization of rriT over unconstrained trial LSP momentum p^. Differentiating m"^ by 
p^, one finds 

This implies that rriT has a stationary value when the trial LSP momentum satisfies 

= = ^P^'. (14) 

This stationary point actually corresponds to the global minimum of tjit for given values 
of m^is and m^: 

(m^) . =771^1 x, = myis + my, (15) 

which is called the unconstrained minimum of the transverse mass 1141. 



niT 




Trial LSP momenta 



Figure 2: A balanced solution. 
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Trial LSP momenta 



Figure 3: An unbalanced solution. 

With the above observation, one can consider two different possibilities for how m^^ 
{i = 1,2) depend on trial LSP momenta, which are depicted in Fig. [2] and Fig. [31 The 
first possibility depicted in Fig. [2] is that m^'' > miji^ for both i when the trial LSP 
transverse momenta take the value giving the unconstrained minimum of m^ji^ {j ^ z), i.e. 

(1)1 ^ (2)1 (2) , 

L'^(l)-_„''"(i)_„''"(2)_-x(2) > rrirp x(2)_-x(2) — "^OTs+"^X; 

I yJrp yJj, \Jrp Pj. I pj, — Pj. 



where p^^'^ is the trial LSP transverse momentum giving the unconstrained minimum of 

(«) 

~x(j) / vis(i) I (i) 

In such case, the minimum of maxjm^'*, m^''} over possible trial LSP momenta is given 
by a balanced mT2 solution [T3] 

771^^2^= min [mJ^M, (17) 

.X(i) 



where the minimization is performed over pj, satisfying 

„X(1) I „X(2) _ visil) vis{2) /.o>, 

Pt + Pr - ~Pt ~ Pr • 



m5!^(p^''^^\p^^^\mj,2,m^) = mP {p^^"^^\p^^^\mP^,m^), 
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The condition ( fT6l) might not be satisfied as in the case of Fig. [3] for which 

'^t''L^(1)- ^^is(2) -x(2) < '"^T'*L>^(2)-f,x(2) = + (19) 

'Pt ~ "r "t "t "t ~Pt 

where p^(^) is the trial LSP transverse momentum giving the unconstrained minimum 
of rrirp . In such case, 112x2 is obviously given by the unconstrained minimum of ruj, . 



Solutions obtained in this way are called the unbalanced 112x2 solution [13] , and given by 

m^f'^' = = m£ + (z = 1 or 2). (20) 

Note that unbalanced solution can accompany a balanced solution as Fig. [3] for which 
the crossing point of m^^ and m^'* corresponds to a balanced 712x2 solution according to 
the definition (|T7|) . However obviously such balanced solution can not be a genuine mx2 
defined as ([5]). 

Recently, analytic expression of the balanced 1^x2 solution for generic event has been 
derived in Ref. [15], which we briefiy recapitulate in the following. For this, let us rewrite 
^^20 as 



T2 



min [max{(ai +/3i)2,(a2 + /52)^}] , (21) 



/3i+/32=v^A-{ai+a2) 



(3l=Pl=ml 



where = {V'^Y — V-Vis the scalar product of the (1 + 2)-dimensional momentum 
vector = {V^, V), and the (1 + 2)-dimensional momenta (y'^,l3^ are given by 

ai — [rjj, ), ^2 — [-t^x 'Pt J; 

/3f = iE^^'\pf\ = iEf\ )). (22) 

These momenta are related as 

< + /3f + < + /?2^ = v^A'^ = v^ (1,0,0), (23) 

where ^/s corresponds to the total transverse energy of the event. The square of the 
balanced mx2 solution corresponds to the minimum of (ai + Pi)^ over all possible values 
of ^/s, subject to the following constraint: 

Pi = Pi = ml, (a, + P,y = {a2 + P2Y. (24) 

The procedure to obtain the balanced 7717-2 solution is first to solve for /3f and P2, 
and next minimize the resulting (ai + Pi)"^ over all allowed values of y/s. Then, one finds 



"'■Throughout this paper, we ignore the initial state radiation, so the total missing transverse momen- 

j- • • u miss I vis(l) , vis{2)'. 

turn IS given by p ji = ^(Pt + Pt )• 
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(r^'lY = ml + AT 



+ 



\ 



( 



\ 



4m^ 



2At - I m. 



(2) 



(25) 



where At is the Euchdean inner product of the two (l+2)-dimensional visible momenta 
ai and ai^: 



A. 



a°a2 + tti ■ 0*2 

r-ivis(l) r-ivis(2) , visil) vis(2) 



(26) 



and we have rewritten the resuh of [15] in a form convenient for our later discussion. 

The mT2 solution obtained above has invariance properties which will be useful for 
the derivation of the possible range of 777,^2 • First of all, the transverse masses m^^ 
and thus 771^2 are invariant under arbitrary independent longitudinal Lorentz boost of each 
mother particle in the pair, which is obviously true as both Ex and p^ are invariant under 
longitudinal boost. Furthermore, 771x2 is invariant also under the following transformation 



of the (1 + 2)-dimensional momenta af = {E, 



vis{i) 
T 



vis(i) \ 

Pt \ 



and 



"2 



At 



-V )a: 



(27) 



where A(^('u) denotes the (l+2)-dimensional Lorentz transformation matrix for 2-dimensional 
boost parameter v. The condition (1161) and the relation (1181) are covariant under the above 
transformation, i.e. if (fT6l) or (fT8|) is satisfied by some (1 + 2) -dimensional momenta, it is 
satisfied also by the transformed momenta. (In appendix A, we provide a more detailed 
discussion on the covariance of f|T6|) and f|T9|) under the transformation fl27|) .) As the 
transverse mass itself is obviously invariant, this means that the balanced solution m^g' 
defined as f[T7|) is invariant under the transformation flTTI) . which can be confirmed also 
by the explicit solution (125|) . As for the unbalanced solution, (fT4|) and (fT9!) are covariant, 
so the unbalanced solution m^^^ defined as (!20l) is invariant also. 

In fact, for generic p"'**^*), the transformation (1271) of the (1 + 2) -dimensional vector 
{Ex, Pt) does not have a direct connection to a true Lorentz transformation of the 4- 
momentum (£", p). However, if both p^**(^) and p*'*'^^^) are in transverse direction, it can 
be identified as a true Lorentz transformation corresponding to a back-to-back boost of 
the mother particle pair in transverse direction. We thus conclude that mx2 for visible 
momenta satisfying invariant under a back-to-back boost of 

the mother particle pair in the direction along the transverse plane T. 



8 



Let us now examine the possible range of mT2, in particular its maximal value. Our 
starting point is the following theorem which will be proved in appendix A: mT2 of any 
event induced by mother particle pair having a vanishing total transverse momentum in 
the laboratory frame is bounded above by another of an event induced by mother 
particle pair at rest. More explicitly, for generic p''**^*) measured in the laboratory frame, 

mT2{pT^'\m^:Lm^) < m^,2(q"^«, m« , mj, (28) 

where q''**(*) is the Lorentz boost of p"*(*) to the rest frame of the z-th mother particle, T' 

is the plane spanned by q''**^^) and q^**(^), so q''**^*) = q^'/^*-* by definition, and the equality 

in the above bound holds when T = T'. Note that T is a fixed transverse plane which 

is independent of the event momenta, while T' varies following the rest frame momenta 
q«s(i) 

By definition, q™*(*) corresponds to the total visible momentum for the decay of the 
i-th mother particle, $j — visibles + x?) measured in the rest frame of and thus its 
magnitude is uniquely fixed by the mother particle mass m, the invariant mass m^*^ of 
the visible part, and the true LSP mass fri^o: 





1 - 


q«s(i) 


2m - 



m + m, - m-o m - m:,-' - m-o 



VIS/ Xi / \ ^ VIS J Xi 



1/2 

(29) 



As a consequence, mT'2 can be described by the three event variables, rn^^is^rn^yis and ^, 
where 9 is the angle between q"*(^) and q"**(2). 

m^,2(q^^^«,mi,mj ^ ^(mi,^,m,). (30) 

Combined with (!28|) . this leads to 

vis(i) (j) 



m?ff (m, 



r / visH) ii) \"| 

max [mT2(Pr ,K/.'"^x)J 



{all events in the lab frame} 
< max [^(mit^,mj] =^-^^(mJ. (31) 

It is always a possible event that both q"***^^) and q'"*^^) are along the direction of the 
transverse plane T, i.e. T' = T. For such event, the invariance of mT2 under longitudinal 
and back-to-back transverse boosts implies 

/ vis(i) (i) \ I vis(i) (i) \ /ooN 

mT2(qT \m'Js^ m^) = mT2(Pr , ^x) (32) 

where now p"*(*) are the momenta obtained by arbitrary back-to-back transverse boost 
along T = T' and subsequent arbitrary longitudinal boost starting from q''**(*). For any 
q"«(«)^ one can choose an appropriate boost for which p''*'^(*) correspond to the visible 
momenta of some event observed in the laboratory frame. This means that for any q''**(')^ 
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there exist some laboratory events whose mT2 is same as the corresponding jF(m|,*g, 6, m^), 
therefore the bound fl3T|) actually corresponds to an identity as 



(33) 



In appendix B, we will show 



< for any m|,*^g, m^, and < ^ < vr, 



dm. 



0=0 



< for < m^o and any m^j 
> for rriy > m^o and any m. 



(i) 



(34) 



and thus the global maximum of JF over the full 3-dimensional event space of {m^^^^, 6} is 
given by 



(35) 



where 



VIS 



"^OTS 5 '"-VIS 



(2) 

ms 

(2) 



- 0, TTiv) 



(36) 



with 



(37) 



In this section, for simplicity we limit ourselves to the one-dimensional event space with 



m 



(1) 



ml/g and ^ = 0, while leaving the discussion of full event space in appendix B. 



For m 



(1) 



(2) 



the unconstrained minima of m^^ have the same value, so mj'2 is 



always obtained as a balanced solution. One can then use the balanced solution fl25|) for 



JF(m, 



(1) 



(2) 



2m 



+ 



= m^is, 9 = 0, m^) 

[m^ - mli^ + rn|o)^ + Afh?{rn?^ - m|o) 



1/2 



2m 



where m is the mother particle mass, m^o is the true LSP mass, and we have used (l29l) for 

m as required. 



|q««(«)|_ Note that ^{m„i.^,m 



m, and thus 



The function ^{myis,m^) has an interesting feature. From 



m^ = m^o^ 



dm„i, 



m 



Xi 



)2p _|_ 477^2(^7722 



Xi 



(39) 
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one easily finds 



dJ^ r < if < m-j^o 



9m„i., I > if > m^o. 



(40) 



This corresponds to the special limit of the general result (134|) . and yields 



^ ^ J^^f^^(m^) = J^(m^i, = m^f ,m^) if > m^o. ^ ^ 



In the above, we have obtained m™2^ from the balanced solution for m\^\. = m\ 



(1) _ (2) 

However, this does not mean that m™!^ is obtained only by the balanced solution. As we 
will see in the next section, m™|^ for > m^o can be obtained also by unbalanced mT2 
solution. In some of such case, the balanced solution giving m™2^ has p™*^ = 0, and thus 
is eliminated by the event selection imposing a nonzero lower bound on |p^*^^| when one 
constructs m™2^ from collider data. On the other hand, the unbalanced solution giving 
m™2^ has a large p™**, so plays an crucial role for the construction of m™2^ from collider 
data. 

If the decay product of each mother particle contains only one visible particle ^^, then 
rriyis is fixed to be m^, and m™|^ is given by 

,2 



(m^ — + ^"^o) + 4m2 (m^ — ""^lo) 



(42) 

2m 2m 
Usually ip is much lighter than the LSP, so one can take the approximation m^ ^ 0, for 
which 



~ 2 2 

m — mtg 



^T2 y^x) - ^ 



im? — m^o 



(43) 



A simple example giving this form of m^g^ would be the squark pair decaj{§, qq — > Q'XiQ'Xi) 
whose will be discussed in more detail in the next section. 

If the decay product of each mother particle contains more than one visible particles, 
m™|^ has an interesting feature which was first noticed in [6]. In such case, m^j^ can vary 
from m™™ to m™^'^, and so m™|'^(m^) at m^^ > m^o takes a different functional form from 
the one at m^ < m^o. As a consequence, m™!'^ has a kink structure, i.e. a continuous but 
not differentiable cusp, at m^ = m-j^o: 

(c;^rV^"^x)»^x="^,o ^ ^ (mg^r)^ - {rnJ^? ^ ^ 

{dJ^^-^/dm^)m^=m.o m2 + ^2 _ (^max)2 ' ^ > 

-^1 Xl 



^Throughout this paper, we do not distinguish squark from anti-squark as wc are considering only the 
kinematics. 
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which becomes sharper when m™f^ becomes larger for a given value of m™lf. 

If visible particles in the decay product are much lighter than LSP, one can set m. 
0, and then 



~ 2 2 

^ + 

2m 



— m?n 
xV 



2m 



+ m2 . (45) 



On the other hand, even for a given number of visible particles in the final state, m^2^ 
can have a different value depending upon the intermediate stage of the decay process. If 
the decay process visibles + does not involve any intermediate on-shell particle 

lighter than the mother particle ^j, one finds 

'^«r = ^-"^x?' (46) 

which results in 

^--(m J = + (m - m^o) . (47) 

However, if visible particles are produced by a chain of decay processes involving inter- 
mediate on-shell particle(s), m™^^ has a smaller value. For instance, if the decay chain 
involves one intermediate on-shell particle with < m, e.g. $.j — + visible 
+ more visibles, the maximal value of m^*^^ for the final state is given by 

<r = ^ (48) 



for which 



-^r^(-x) = r^(i-^) + ^(i- ' 




m^ 2 m^ 



+ 



i 





-V -T^-^ +"^x- 49 
m^ 2 m^ ' 



In Fig. m we depict the behavior of m™2^(m^) for the case without any intermediate 
on-shell particle {m^ > m) and also the case of decay chain involving an intermediate on- 
shell particle {m^ < rh). As the value of m^f for the case without intermediate on-shell 
particle is bigger than the value of m™^^ for the other case, the kink structure of Fig. |4](a) 
is sharper than that of Fig. 111(b) as anticipated in (jS]). A simple example of the process 
that each mother particle is producing more than one visible particles is the gluino pair 
decay: gg QQXiQQXiy whose mT2 will be discussed in more detail in the next section. 
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■"z "z! "z 



Figure 4: m™!"^ for (a) m,^ > m, (b) < m 

3 Features of squark and gluino 771^2 

In this section, we discuss in more detail the 772^2 of two specific processes, the decay of 
pair-produced squarks, qq — > qXilXi: example of the case that each mother particle 
decays to one visible particle and one invisible LSP, and the decay of pair-produced 
gluinos, gg qqXiQQXi^ an example of the next case that the decay product of each 
mother particle contains more than one visible particle. We also perform a Monte Carlo 
LHC simulation for some superparticle spectra to examine how well can and m^g'^ 
be constructed from real collider data. 

3.1 Squark mT2 

Let us consider mT2 for the process in which a pair of squarks are produced in proton- 
proton collision and each squark decays subsequently into one quark and one LSP: 

pp-* qq-^ qXiQXi, (50) 

where q denotes the 1st or 2nd generation quark. A characteristic feature of this process 
is that = '^iil = "'^g; and thus m^"* and m^^ have the same unconstrained minimum: 

("^?)mi„ = ("^?Vin = + ^l- (51) 

Then, 771x2 is always obtained as a balanced solution as shown in Fig. [2l and the resulting 
balanced solution fl25|) is simplified as 

mT2 = ml + AT + MAT-ml + 2ml^(^AT + miy (52) 
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giving 



mT2 



where 



At = E, 



+ 



vis(l) j-,vis{2) 



At — 171'^ + 2m^ 



(53) 



vis{l) vis{l) 



Pt 



T-ims(l) T-i«s(2) 



I visll) 1 1 vis{2) I n 

IPt IIPt I "^ost^ 



(54) 



for E, 



vis(i) 



T 



Pt 



ms(j) 12 



(^ = 1,2). 

As was discussed in the previous section, mT2 of any event induced by mother particle 
pair having a vanishing total transverse momentum in the laboratory frame is bounded 
above by another 771x2 of an event induced by mother particle pair at rest. For a squark 
pair at rest, the magnitude of quark momentum from each squark decay is given by 



visit) I 



2mn 



((mg + 77lqf - m|o) {{77lq - 77lqY - m|o) 



1/2 



(55) 



It is then straightforward to find that 771t2 has a maximum at ^ = for visible momenta 
in transverse direction, giving 



771t2 



m? - m|o + ml 

277lq 



+ 



\ 



[{77lq + 77lqY - m|o) {{niq - 77lqY - m?o) 



4m? 



771X 



(56) 



as obtained in (l42l) . In the limit when the quark mass is negligible, this expression of 
m™|^ is further simplified as 



T2 



771] - m?o 



2m^ 



2 2 \ 2 



(57) 



As the above m™2^ has been obtained from a momentum configuration in which the 
two quarks produced by a squark pair at rest are moving in the same direction, one might 
worry that such two quarks might not be identified as two separate jets in real collider 
data with realistic jet reconstruction. However, the quark momenta from the back-to- 
back transverse boosted squark pair are generically not in the same direction, and thus 
the boosted events can provide two separate jets, while giving the same value of m™2^. 

To see explicitly some features of the squark 771x2, we have performed a Monte Carlo 
analysis for a SUSY parameter point in mirage mediation scenario [16], which gives 



77iq = 697 GeV, 



771^0 
Xi 



344 GeV. 



(58) 
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A Monte Carlo event sample for the signal pp ^ qq ^ QXiQXi has been generated in 
partonic-level using the PYTHIA event generator [17]. The mT2 values for the event 
sample were then calculated with a numerical code [TS] implementing the minimization 
over the trial LSP momenta. Fig. shows the resulting distribution for a trial 
LSP mass = 10 GeV. The distribution shows a sharp edge at 771^2 = 527.8 GeV, which 
is very close to m™|^ = 527.4 GeV obtained from the analytic formula (!57l) . Finally, m™|'' 
as a function of the trial LSP mass is shown in Fig. Mjo). Here, the blue curve 
represents the analytic formula (1571) . while the black dots are obtained from the Monte 
Carlo data, which fit very well the analytic curve. 




Figure 5: (a) mT2 distribution with = 10 GeV for the mirage mediation parameter 
point fl58l) . (b) Resulting m™2^ as a function of the trial LSP mass m^. 

Fig. [6] shows the distribution of the opening angle 6 between the two visible quark 
momenta for the events near (a) the lower edge and (b) the upper edge of Fig. [5t^a). As 
anticipated, the events near the lower edge and the upper edge give a distribution with a 
peak at cos 6' = —1 and cos^ = 1, respectively. The rather broad peak in Fig. MJo) is due 
to the non-zero transverse momentum of the each squark, which makes the two quarks 
from the squark-pair decay less aligned in general. 

3.2 Gluino 111x2 

Let us now examine a more complicate process in which each of the pair-produced mother 
particles decays into one invisible LSP and more than one visible particles, for which 
m^2^{m^) shows a kink structure at = m^o as was noticed first in [0] and discussed in 
the previous section in a generic context. As a specific example, we consider the process 
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Cos6 



Cos9 



Figure 6: Event distribution along cos 9 for the events near (a) the lower edge and (b) the 
upper edge region of the Fig. ^a). 

in which a pair of gluinos are produced in proton-proton collision and each gluino decays 
into two quarks and one LSP: 

PP^gg^ QQXiQQXi, (59) 

where again q denotes the 1st or 2nd generation quark. Depending upon whether squarks 
are heavier or lighter than gluino, the gluino decay g QQXi occurs through a three- 
body decay induced by an exchange of off-shell squark or two body cascade decay with 
intermediate on-shell squark. As noticed in the previous section, these two cases have a 
different value of m^f^ = m™^^, thereby the resulting ■m'^2^{m^) has a different functional 
form in the range > m^o. 

For the process ( l59l) . the unconstrained minima of the transverse masses of the decay 

product of each gluino, i.e. {''^t^),^^^ = ^vL + "^x' generically different from each 
other, thereby both of the balanced and unbalanced solutions can appear. 

If the squark mass rriq is heavier than the gluino mass rrig, gluino will undergo the three 
body decay g —>■ qqx\ through an exchange of off-shell squark. In this case, for = 0, 
the total invariant mas^ of the visible part of each gluino decay is in the following range: 

< mi2, < rUg- m^o. (60) 

As discussed in section 2, in order to obtain the maximum of gluino mT2, we can limit 
ourselves to the situation that the two gluinos are produced at rest and all decay products 

^Onc might consider the total transverse mass of the visible part which has the same range. 



16 



are moving on the transverse plane. In such case, the transverse momentum and energy 
of the visible part are given by 



I vis(i) I 

IPt I 



^J {(iT^g + m^ly - ^o) {{m-g - m^XY - ^o) 



E. 



vis{i) 



ml - ml, + (m2j2 



2mg 



(^ = 1,2), 



2mn 



(61) 



Then, the balanced mT2 solution ( !25l) is unambiguously fixed by the four variables: m^2) 
m[2, 6 = the angle between p^****-^^ and p^^^^'^\ and the trial LSP mass m^. 

For given values of m|,*g and m^, the maximum of the balanced mT2 occurs at ^ = 
(see the appendix B for a detailed proof.) Since we are mainly interested in m™|^, we will 
focus on the configuration with ^ = in the following. From 0251) . we obtain 



m^2\e=o = "^x + ^ 



Ami 



2A- m. 



m 



(2) 



(1) (2) A' 



A? [ Tn^^^m^^^ J I , 



(62) 



where 



A = A\, 



j-ivis{l) j-,vis{2) I vis{l) 1 1 vis(2) i 

rVrp tjrp + I Y>rp \ \ ^irp | , 



(63) 



with |p^**^*^| and E™^^"^"^ given by (l6Ti) . If we further restrict the momentum configurations 
to satisfy 



(1) (2) _ 



(64) 



the balanced mT2 solution is simplified as 



mT2\e=Q 



2mh 



+ 



{[vfig + m^i.Y - m\) {{m-g - m^i^y - m^o) 



4m? 



— + m2. 



(65) 



Then, from (HH), we obtain 



m??r(m,) 



m??r(m, 



2mg 



\ 



2 2 \ 2 

^ ' + if < m^o , 



2mg 



mg — m^o ] + m^. 



if m^ > m^o. 



(66) 
(67) 
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Figure 7: Extreme momentum configuration providing m™2^ as a balanced solution for 
(a) < m^o and (b) > m^o. 

Fig. Wi^) shows a momentum configuration providing the m™2^ of fl66|) . In this config- 
uration, two gluinos are produced at rest, and each gluino subsequently decays into two 
quarks moving in the same direction (i.e. m^2 = ''^iis — 0) LSP moving in the 

opposite direction. Furthermore, two sets of gluino decay products are parallel to each 
other (i.e. ^ = 0) and all of them are on the transverse plane with respect to the proton 
beam direction. This configuration is the first example of extreme momentum configu- 
ration considered in Ref. [6]. Although this corresponds to the simplest configuration 
providing the m™|^ of (l66!l . it might not be useful for constructing m™2^ from real collider 
data as all quarks are moving in the same direction, so that they can not be identified as 
separate particles due to the finite jet resolution. 

This difficulty of jet resolution can be partly avoided by the back-to-back transverse 
boost of the above extreme configuration, giving the same value of m™2^(m^). In the 
back-to-back boosted configurations, the two di-quark systems are not moving in the 
same direction in general, so that can be distinguished from each other in real collider 
event. However, the two quarks in each di-quark system are still aligned to each other. 
In real collider data analysis, two aligned quarks cannot be identified as separate jets 
with realistic jet reconstruction, which will eliminate the events which involve the quarks 
moving in the same direction. As the true maximum of comes from such momentum 
configuration, any realistic jet reconstruction will cause a systematic shift of m™2^ to 
a lower value when one tries to construct m™|^ from real collider data. Our analytic 
expression fl25l) for the balanced solution provides information on how sensitive 171^2 
is to the angular separation of the involved quarks, with which one can estimate the 
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uncertainty of m™|^ caused by the jet resolution cut: 

-^[1-^]{ArY <Oil)%, (68) 

where AR = y^A0^~+~Ar/^ ~ 0.5 represents a separation of two quarks in azimuthal angle 
and pseudorapidity plane. This indicates that the systematic shift of m™|^ due to the 
finite jet resolution is negligible, which we have confirmed by an explicit Monte Carlo 
analysis. 

The momentum configuration of Fig. [7](b) provides the m^2^ of (1671) . which was 
considered in [B] as the second example of extreme momentum configuration. Here, gluinos 
are pair produced at rest, the two quarks from each gluino are back to back to each other 
("^i2 = "^Ss = ^9 ~ "^Xi')' while the LSP is at rest. In this case, the angle 6 is not well 
defined because p^**^^'' = p^***-^-* = 0. Also p^**'^ = which is true even after a back-to-back 
boost of the system. Such momentum configuration will be useless when one constructs 
m™2^ from collider data as one typically uses an event selection cut imposing a lower 
bound on |p™''^^|. However, as we will see shortly, there exist momentum configurations 
yielding fl67|) while having a sizable |p™^*|, so that the m™|^ of fl67j) can be constructed 
from collider data even under a proper cut on |p™'^^|. 



''T2 




Figure 8: Extreme momentum configuration providing m™|^ as an unbalanced solution 
for rriy. > m^o. 

A momentum configuration, which provides the m^g^ of fl67|) with a sizable |p??*^'^|, 
is shown in FiglHl In this configuration, two gluinos are produced at rest. The first 
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gluino produces a di-quark system with m|,2 = and one LSP, while the second gluino 
produces a back-to-back di-quark system with m„j^ = rrig — m^o and one LSP at rest. 
For the second gluino decay set, the visible momentum p^****-^'* = 0. Thus, unconstrained 
minimum of the second gluino transverse mass, = {nig — m^o) + m^, occurs when 
trial LSP momentum p^^^'' = 0. On the other hand, the first gluino decay product has 
|p«s(i)| _ _ ^2^j)/2mg and p?.*'^^ = — p^*'^'-^^ for pJ.*'^'' = 0. Then, the corresponding 
transverse mass of the first gluino decay is given by 



(1) - "^1? 

= 

^ 2m~g 



\ 



2 2 \ 2 

2mg 



+ m2 . (69) 



Therefore, > m^'' if > m^o, though m^'' is at the unconstrained minimum value, 
so that we have unbalanced mT2 solution, i.e. mT2 = ("^g ~''^x?) "'"'^x ^^^^ momentum 
configuration. The same unbalanced mT2 solution is obtained for the momentum config- 
uration with other values of m|^2 because those cases also give m^^ > rrSp when m!p 
is at the unconstrained minimum. Such momentum configurations and the back-to-back 
transverse boosted ones would have a sizable |p™'^'*|, so can be used to determine m™|^ 
from real collider data. 

To summarize the extremal features of mj^2 for the decay of gluino pair when > m^, 
the maximum of over all events is given by 



m? - rriio 



2mg 



\ 



2 2 \ 2 

+ rnl if niy < m^o, 

2mg I X xi' 



rag — rri^o^ + if > m^o (70) 



as obtained in (145!) and (1471) in more generic context. Thus, there is a level crossing of 
m™|^ at m-^ = fn^o, yielding a kink structure as shown in Fig. |4](a). If such m™|^-curve 
can be constructed from collider data, which will be examined in the next subsection, this 
kink structure will enable us to determine the true LSP mass m^o and the gluino mass 

Xi " 

nig = m™|^(m^ = m^o) simultaneously. 

To see explicitly the extremal features of gluino niT2, a Monte Carlo event sample 
of the signal pp ^ gg QQXiQQXi has been generated in partonic-level, for a SUSY 
parameter point in a minimal anomaly mediated SUSY-breaking (AMSB) scenario [19], 
which gives 

m-g = 780 GeV, m^o = 98 GeV, (71) 

with a few TeV sfermion masses. The niT2 values for the event sample were then calcu- 
lated. Figj9] (a) and (b) show the resulting niT2 distributions for trial LSP mass = 10 
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Figure 9: 'mT2 distribution with (a) = 10 GeV and (b) = 350 GeV for the AMSB 
parameter point (17T|) . 




Figure 10: m™2^ as a function of the trial LSP mass for the AMSB parameter point 

GeV and 350 GeV, respectively. On the figures, hatched histogram corresponds to the 
balanced values, while black histogram to the unbalanced ones. As anticipated, one 
can notice that for = 10 GeV, which is smaller than true LSP mass, the endpoint of 
the distribution is determined by the balanced mj'2 solutions, while both balanced 
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and unbalanced mT2 solution contribute to the endpoint region for = 350 GeV, which 
is larger than true LSP mass. Finally, m™|^ as a function of the trial LSP mass is 
shown in Fig. [TDl Here, the blue and red curves represent the analytic formula (1701) . while 
the black dots are obtained from the Monte Carlo data, which fit very well the analytic 
curves. 

Let us now consider the case of lighter squark, rriq < rrig, for which the following 
cascade decay is open: 



qq 



(72) 



where the squark in the second stage is on mass-shell. The main difference between this 
two body cascade decay and the three body decay is that the total invariant mass of the 
visible part takes the range: 



n ^ (1) (2) 
< mV^, m^^ 



< 



m?)(m? 



mi 



(73) 



As m'^f^ has a smaller value than the three body decay case, while is same, the kink 



structure is weakened as was anticipated in (jH]). The maximum of mj'2 for m^^ < m^o 
takes the same form as the case of heavier squarks, while it is changed to a different form 
for m^ > m^o: 



m 



T2 



m??r(mj 



+ 



mi 



2mg 



\ 



mt — m - 



2mn 



2 mi 2 



' + ml if m^<m^o, 



Xi 




mi 



+ m| if > m^o (74) 



as obtained in (149|) . The m™2^-curve for lighter squarks is depicted in Fig. ID^b), which 
shows again a kink structure at = m^o, although milder than the case of heavier 
squarks. Note that the m^f^-curve for the range m^ > m^o depends on the squark mass 
also, so it can determine the gluino mass, the LSP mass and the squark mass altogether. 



3.3 Construction of gluino mT2 from collider data 

In order to check the experimental feasibility of measuring superparticle masses using 
the kink structure of the gluino m™2^, we have generated Monte Carlo event samples of 
the SUSY signals at LHC by PYTHIA [T7] for several SUSY breaking schemes yielding 
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different patterns of superparticle spectra. We have also generated the SM backgrounds 
such as tt, W/Z + jet, WW/WZ/ZZ and QCD events, with less equivalent luminosity, 
in five logarithmic pt bins for 50 GeV < pt < 4000 GeV. The SM backgrounds have 
been also generated by PYTHIA. The generated events have been further processed with 
a modified version of the fast detector simulation program PGS [20], which approximates 
an ATLAS or CMS-like detector with reasonable efficiencies and fake rates. 

For each event, the four leading jets are used to calculate the gluino mT2- For a 
convenience of numerical analysis, we considered mT2 defined in terms of the transverse 
visible mass m^^^^\ rather than in terms of the invariant visible mass m^j^^^, which gives 
the same value of m™2^ as remarked in section 2. Note that m^'^^^^ = rrS^J^ for the extreme 
momentum configurations giving the maximal value of mT2- The four jets are divided in 
two groups of dijets as follows [21]. The highest momentum jet and the other jet which 
has the largest \pjet \ with respect to the leading jet are chosen as two 'seed' jets for 
the division. Here, pj^t is the jet momentum and Ai? = ^/AcjP^^T^Tf denotes the jet 
separation in azimuthal angle and pseudorapidity plane. Each of the remaining two jets 
is associated to the seed jet making a smaller opening angle. Then, each of the jet pairs 
constructed in this way is considered to be originating from the same mother particle 
(gluino). If this procedure fails to choose two groups of jet pairs, we discarded the event. 

Because the functional form of the gluino m™2^(m^) depends upon whether the 1st 
and 2nd generations of squarks are heavier than gluino or not, we consider those two 
cases separately. Let us first consider the case of heavier squarks. Superparticle spectrum 
with rriq > rrig can arise from various SUSY breaking schemes, for which the gluino mT2 
takes the form of Fig. Hl^a). For simplicity, here we consider only the case that the 
1st and 2nd generations of squarks are significantly heavier than 1 TeV, e.g. rriq ~ 4 
TeV, so that squarks are not copiously produced at LHCf], while gluino is light enough 
to be copiously produced, e.g. nig < 1 TeV. As an specific example of such superparticle 
spectrum, here we consider a parameter point of anomaly mediation scenario (AMSB) 
with heavy sfermion masses, in which the gluino, LSP and (the 1st and 2nd generation) 
squark masses are given by**! 

AMSB with heavy sfermions : rrig = 780 GeV, m^o = 98 GeV, rriq = 4 TeV. (75) 

For the AMSB point, the production cross section of gluino pair cr{gg) ~ l.lpb. The 
branching ratios of gluino decay are as follows: B{g XiQl) ~ 32%, B{g — > X2II) ~ 3%, 
B{g xfqq') ~ 64%. Thus, gluino mostly decays into lighter chargino (or LSP) plus two 

"if squarks are heavier than gluino, but still light enough to be copiously produced at LHC, the gluino 
mT2 is not a proper observable to determine the gluino and LSP masses since the gluino-pair events 
99 ~^ QQXillXi s-re severely screened by the squark pair events qq qgqg —* qqqxIqqqXi- In such 
case, one can construct the squark mT2 for the squark pair events with six quarks, which would show 
a behavior similar to the gluino m^2^ in the case of lighter squark, from which the squark, gluino, and 
LSP masses can be determined altogether. 

**We also assume tan/3 = 10 and the Higgsino mass parameter ^ > 0. 
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quarks. Being wino-like, the LSP and the hghter chargino are almost degenerate in mass. 
The chargino decay xf ~^ Xi^^^ produces very soft leptons, which cannot be detected at 
LHC In this circumstance, both gluino decays g — * XiQQ.' ^i-nd g Xiqq can be considered 
as 'signals' we are looking for, and the contamination from the small number oi g x^qq 
decay is expected not to be significant. In this work, we assume integrated luminosity of 
300 fh-^ for the AMSB point. 
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Figure 11: (a) Gluino distribution with = 90 GeV for AMSB with heavy 
sfermions, and (b) m™2^ as a function of for AMSB with heavy sfermions. 

To obtain a clean signal sample for the gluino mT2, we have imposed the following 
event selection cuts on the SUSY and SM event samples. 

1. At least 4 jets with Pti,2,3,4 > 200, 150, 100, 50 GeV. 

2. Missing transverse energy EJp^'"^^ > 250 GeV. 

3. Transverse sphericity St > 0.25. 

4. No b-jets and no leptons. 

Using the event set passing these selection cuts, we calculate the gluino 771^2 for various 
values of the trial LSP mass m^. Fig. [TT] (a) shows the resulting gluino mT2 distributions 
for the AMSB with = 90 GeV. Fitting the distribution with a linear function with a 
linear background, we get the endpoint value 

AMSB : m??f (m^ = 90) = 778.2 ± 2.2 GeV (76) 

The edge values of obtained in this way are shown in Fig. [TT] (b). Blue and red lines 
denote the theoretical curves obtained from fl70|) . Fitting the data points to these curves. 
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we obtain the following gluino and LSP masses: 

AMSB : m-g = 776.5 ± 1.0, m^o = 94.9 ± 1.4 GeV, (77) 

which are quite close to the true values in (1751) . This demonstrates that the gluino mr2 
can be very useful for measuring the gluino and the LSP masses experimentally in heavier 
squark scenario. 

Let us now consider the case of lighter squarks, rUq < rUg, for which the cascade decay 
g ^ qq qqXi is open. As an example of superparticle spectra with lighter squarks, we 
choose a parameter point (SPSla [22]) of mSUGRA schemes, which provides 

mSUGRA with light squarks : nig = 613, rriq = 525, m^o = 99 GeV. (78) 

For this mSUGRA point, the production cross sections for gg, gq and qq pairs are <j{gg) ~ 
4.2 pb, <y{gq) ~ 21 pb, and cr(gg) ~ 9 pb, respectively. The branching ratio of the signal 
decay chain, i.e, g ^ qq Xiqq is B{g Xiqq) ~ 40%, while corresponding branching 
ratios to X2, and xf are B{g — >• X2qq) ~ 7%, and B(g xfqq') ~ 14%, respectively. 
Here, we assume 30 fb~^ of integrated luminosity for the mSUGRA point. 

Similarly to the above AMSB case, we have imposed following event selection cuts: 

1. Missing transverse energy EJp^''^'^ > 250 GeV. 

2. Transverse sphericity St > 0.25. 

3. No b-jets and no leptons. 
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Again using the Monte Carlo events passing the selection cuts, we evaluated gluino mT2- 
FigHS] (a) shows the mT2 distribution for trial LSP mass = 90 in mSUGRA scenario 
with light squark. Even with the above selection cuts, the events from the gq pair pro- 
duction largely contribute to the 171^2 distribution. The contribution from gq pair events 
provides rather similar shape of 771^2 distribution to the one from gg events, but the 
maximum of mT2 from gq events is still smaller than the one from gg events. 
Fitting the edges of these distributions, we find the endpoint values: 

mSUGRA : m?!f (m^ = 90) = 610.8 ± 2.1 GeV. (79) 

Figll2] (b) shows m™|^ as a function of m^. Fitting the data points to the analytic 
expression ([74]), we obtain 

mSUGRA : = 611.7 ±2.8, = 519.9 ± 2.8, m^o = 96.3 ± 8.1 GeV. (80) 

which are again quite close to the true mass values in (ITHI) . We have performed the same 
analysis for other superparticle spectra with rrig < rrig, e.g. a parameter point of mirage 
mediation scenario [16], and found that the gluino mass, squark mass and LSP mass can 
be determined with a similar accuracy. 

Here, we emphasize that the above results include only the statistical uncertainties. 
There should be various systematic uncertainties associated with the choice of fit function 
and the fit range to determine the endpoint of distribution, which would affect the 
result. Study of such systematic uncertainties, however, is beyond the scope of this work. 

4 Conclusion 

In this paper, we provided a detailed study of the collider observable mT2 applied for 
pair-produced superparticles decaying to visible particles and a pair of invisible LSPs. In 
particular, we have derived the analytic expression of the maximum of mT2 over all events 
(m™!""). It is noticed that if the decay product of each superparticle involves more than 
one visible particle, m™2^ being a function of the trial LSP mass has a kink structure, 
i.e. a continuous but not differentiable cusp, at = true LSP mass, which can be 
used to determine the mother superparticle mass and the LSP mass simultaneously. The 
sharpness of the kink structure depends on whether the full decay process involves an 
intermediate on-shell particle (lighter than the mother particle) or not. In case without 
any intermediate on-shell particle, the kink structure is sharper. In other case with an 
intermediate on-shell particle, although the kink structure is weakened, the m^l'^-curve 
can be used to determine the intermediate particle mass also. 

We also performed a Monte-Carlo study of the gluino mT2 for some superparticle 
spectra in order to examine how well m™2^ can be constructed from collider data. The 
result of our study indicates that the kink structure of m™!"" can be quite useful for the 
determination of superparticle masses in many cases, and determine the mother particle 
mass and the LSP mass quite accurately in some cases. 
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Appendix A 

In this appendix, we show that 771^2 of any event in the laboratory frame is bounded 
above by another 771x2 of an event induced by mother particle pair at rest. Let us consider 
generic event induced by a symmetric decay of mother particle pair: 

Xi+ visible particle(s), (81) 

and let p"'**^*) (i = 1,2) denote the total visible momentum of the decay product of <l>j 
measured in the laboratory frame. The corresponding is determined by the visible 
transverse momenta p^**^*-* and the visible invariant masses m^*^^ as defined in ([5]): 

7nT2{pT^'\m^Lm^) = w,, [max {mjl^ m^^^}] , (82) 

where 



= ^7711 + (m« )2 + 2i?^^^«E^« - 2p;^^(^) ■ p^«. (83) 

As the first step, let us perform independent longitudinal boost of mother particles to 
make the mother particle pair to move back-to-back in transverse direction. Note that one 
can always make such longitudinal boost for mother particle pair having a vanishing total 
transverse momentum in the laboratory frame, and also the required longitudinal boost 
of $1 is generically different from the one of $2- Let p'™*(*) denote the visible momentum 
after such longitudinal boost of mother particles. As pr and Et = \/ 771^ + | Pt p are 
invariant under longitudinal boost, we obviously have 

mT2(pr^^"^!2,^x) = mT2(p'"^«,m« ,m,). (84) 
To proceed, let us consider the generalized invariant mass m/2 defined as follows: 

m,2(p"^(^\m«,m,)= ,min [max {m«,m(2)}] , (85) 

where m^*^ is the trial invariant mass of the mother particle $i obtained for the trial LSP 
mass m^: 



m(*) = ^m\^ {m^vly + 2E«^(*)Ex(*) - 2p"*(*) ■ px(i), (86) 
and the minimization is performed over the trial LSP momenta satisfying 

pX(l) + pX{2) ^ _pv^s(l) _ p«s(2)_ ^g^) 

As m^"* < m^*-* for arbitrary visible and trial momenta, one immediately finds 

mr2(pr^^ma,m,) < m,2(p'"^(^ m« , mj, (88) 
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where the equahty holds for p'"*(*) in T-direction. Note that m/2 can be considered as an 
(1 + 3)-dimensional analogue of the (1 + 2)-dimensional mT2- 

The global minimum of m^*) over the unconstrained p^(*) is given by 

K^)min = "^!.t + '^X, (89) 

which occurs when 

pX» = E'^i^^visii) i^vis(i) ^ ^^p«^» . (90) 

Like mT2i the generalized invariant mass m/2 is also given by either a balanced solution 
or an unbalanced solution. If m^*) > m'^-'^ for both i when the trial LSP momenta take 
the value giving the unconstrained minimum of m'^-'^ (j ^ i), i.e. if 



^\pX(l)=-pvis(l)_pV^s(2)_px(2) — | pX{2) =px(2) = ''^Li + "^X ' 

|pX{2)=_p«is(l)_pi.»s(2)_pX(l) — ''|pX{l)=pX(l) ^ ■'^12 + "^X' (91) 

where 

the corresponding mn is given by a balanced solution: 

mi2 = m}if= min [m^^)], (92) 

{pX(l)+pX(2)=_pi)is(l)_p«»3{2)^^{l)=m(2)} L J 

where the minimization is performed over p^^*) satisfying 

^(i)(p-(i),px(i),^{i),^j = m(2)(p-(2),px(2)^^(2)^^j^ 

pX(l)+pX(2) _ _p«^(l) _p««(2)_ (93) 

On the other hand, if m*^*^ < m'^-'^ for any z when the trial LSP momenta take the 
value giving the unconstrained minimum of m*--'^ (j 7^ 2), i.e. if 

^^^^|px(o=px« < ^^^■1pxu)=px(.) (jVO (94) 
for the trial LSP momenta given by 

pxOO = px{i) = m^p--(i)/m!,2, 

pX{2) _ _p«s{l) _ p«s(2) _ pX(l)^ (95) 

or 

pXO-) = pX(2) ^ ^^p««(2)/^(2), 

pX(l) = _p-{l)_p-{2)_pX(2)^ (gg) 
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the corresponding 17112 is given by an unbalanced solution as 



mj2 



m 



unbal 
12 



m 



(i) 



+ (j = 1 or 2). 



(97) 



In section 2, we have noticed that mT2 is invariant under the back-to-back boost of 
p««(*) along the direction of the transverse plane T, if both p"*(^) and p''*'*^^) are in the 
direction of T. It is in fact straightforward to show that m/2 is invariant under back- 
to-back boost in general direction for general p"*(*), which corresponds to the (1 + 3)- 
dimensional version of the back-to-back boost invariance of mT2- To show this, let us 
first note that the invariant masses m*^*^ and the relations (l9Ti) and (1931) are invariant or 
covariant under the following back-to-back Lorentz boost: 



ai — > 

a^2 ^ K{-v)a''2. (3^~^K{-vW2, 



where 



vis{i) vis{i) 
) P 



(98) 

are the visible 4- momenta, /5f = {E^^'^\ p^^*^) are the trial LSP 



4-momenta, and A^(v) denotes the (1 -|-3)-dimensional Lorentz transformation for generic 
3-dimensional boost parameter v. This assures that m/2 given by a balanced solution 
is invariant under generic back-to-back boost, which can be confirmed by the following 
explicit form of rri^f: 



bal\2 



A 



+ 



v 



4m^ 



2A 



(1) 



m. 



(2) 



(1) (2) 



(99) 



where A is the Euclidean product of the two visible 4-momenta and 02'- 



A = E™^^^'^ _£'"*(2) _|_ p«s(l) . p«s(2)^ (100) 

Similarly, the relations 0941) and (1951) are covariant under the above back-to-back Lorentz 
boost, so m/2 given by an unbalanced solution is invariant also. 

In fact, the covariance of fl9T]) and flM|) under the back-to-back boost fl98l) . which we 
have used to show the invariance of m/2, is not so obvious. The easiest way to see their 
covariance is to consider the boundary between fpT]) and which corresponds to the 
visible momenta satisfying 



^(i)(p-(i)^px{i)^^(i) 



for the trial momenta given by 

,x(i) 



m^) 



m(2)(p-«(2),px{2)^^(2!,m 



(101) 



p 



m^p 



vis{l) 



X(2) 



P 



x(i) 



/ (1) 
pws(l) 



vis{2) 



(102) 
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or by 

pX(l) = _pX(2) _p«s(l) _p«s(2)_ (^Q3) 

If some momenta satisfy fllOll) and (11021) . or fllOip and (11031] , their back-to-back boost 
satisfy also the same relations. This means that the boundary is invariant under the back- 
to-back boost, thus there can not be any crossing of the boundary caused by the back- 
to-back boost. Therefore, if some visible momenta satisfy (19T|) or (!94|) . their back-to-back 
boost satisfy the same condition, so ( jOTl) and (|94l) are covariant under the back-to-back 
boost (l98l) . The same argument can be used to show that (fT6l) and (fT9l) are covariant 
under the (1 -|- 2)-dimensional transformation (1271) . 

In the above, we have argued that mj2 is invariant under generic back-to-back Lorentz 
boost, which leads to 

m,2(p'"^«,mi ,m,) = m,2(q"^«, m^l, m,) (104) 

for q''*^(*) obtained by arbitrary back-to-back boost of By definition, p'"^(*) is the 

i-th visible momentum after the independent longitudinal boosts making the mother par- 
ticle pair to move back-to-back in the direction of T. One can then choose an appropriate 
back-to-back boost for which q*'*''(*) corresponds to the i-th. visible momentum measured 
in the rest frame of its mother particle. 

Let T' denote the transverse plane spanned by q"*(^) and q^*'^^^), and z' denote its 
normal direction. Then, by definition q"^*^*^ = 0, and it is straightforward to minimize 
max jm'-^'*, m*^^-* } over the 2;'-component of the trial LSP momentum: 

min [max jm*'^'', m'-^Hl 

{px(l)+px(2)=_qi.is(l)_qfis(2)} L J J 

min [max{m'^^\m^^H x(i)_„x(2)_ 

[ max \^m^j,} , m^,"* }] , (105) 



mm 



and thus 



m,2(q"^«,m« ,mj = mT,2(q^^^ m« , m,). (106) 

Combining dHl]), ([88]), ffTOD and ffT06D . we finally obtain 

m^2(p"^«,m« ,mj < mr,2(q"^«, m« , (107) 

for arbitrary symmetric decay of mother particle pair having a vanishing total transverse 
momentum in the direction of T, where the equality holds when T' = T. Therefore, 
mT2 of any event induced by mother particle pair having a vanishing total transverse 
momentum in the laboratory frame is bounded above by another mj-2 of an event induced 
by mother particle pair at rest. 
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Appendix B. 



In this appendix, we discuss in detail the 'mj'2 of the events associated with the decay of 
mother particle pair at rest with visible momenta in transverse direction. Such 771^2 corre- 
sponds to mT/2(q*'*'*''*'', m^*g, 771;^), where q™*(*) denotes the z-th visible momenta measured 
in the rest frame of its mother particle, and T' is the transverse plane spanned by q"**(^) 
and q™*(2)_ gin^g |q«s(«)| jg determined as 



|q 



uis(i)| ^ 

2m 



m + m; ) - m-o ] ( [m - m:,/J - m-o 



1/2 

(10^ 



where m is the mother particle mass and m^o is the LSP mass, mj'/2 appears as a function 

of the three event variables m^^^g,m^^^g,6 = the angle between q"*(^) and q"*(^), and also 
the trial LSP mass m^: 

mT'2{q'"'^'\m^:lm^) ^ J-(m«,^,mJ. (109) 

In section 2, we discussed the behavior of over the one-dimensional event space with 
m^2 = 'f^^vis = "^fis and 6 = 0. Here we generalize the analysis to the full 3-dimensional 
event space spanned by {m^*^^,^}, and show 



< for any m^^]^,m^, and < 6 < tt, 



dm:,- 



(i) 
vis 



> for m^ > m^o and any m. 



< for m^^ < m^o and any m^^^^, 



:iio) 



and thus the global maximum of over all {m^-g, 6} is given by 



where 



•^^ = J'i^ = <t,rr^ = <i:,0 = 0,m^, 

^ =^max^^(21 = m-r,^ = 0,m,). (112) 



As discussed in section 2, 771^2 is given by either a balanced solution or unbalanced 
solution, depending upon whether the condition f|T6|) is satisfied or not: 

_ I JF*^^' for q*'*'^^*) in the balanced domain, 
= jrunbai otherwise, 
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where 



for 



bal\2 



+ 



\ 



1 + 



4m^ 



- (^S)^ - irn^' 



r 



unbal 



+ ^vl = 1 or 2) 



:il3) 

:ii4) 
:ii5) 

:ii6) 



with |q«'*(*)| given by ffTOHj) and E'''^''^ = i^|q"^»|2 + (m^J^- Here the balanced domain 
corresponds to the event set satisfying fITBl) . and JF'^'^^ is obtained from fl2S]) with p™**-*-* = 
q««{«)_ Note that q''**^*) = q^'/^*^ and A = A^/ according to the definition of T'. 



B.l T vs. ^ 



Let us show that JF has its maximum at 6' = for given values of m^y]g and rriy. If 



(1) 



(2) 



rrims, one easily finds J-' is given by a balanced solution of the form 



which obviously has its maximum at = for given values of m^/g and m^^. 
To examine the case with m|,2 7^ "^12' consider 



bal\2 



89 



|q«s(i) I |q«s(2) I g-j^ ^ 



where 



A = _£'«''{2) _j_ |q' 



iiis(l) I I vis{2) 



q 



cos6', 



X 



2A 
A^ 



[m. 



W\2 



[m. 



{2)n2 



/ (1) (2)\2 



B 



It is straightforward to find that A > 0, 5 > and C > when m^j^ ^ m 
has an extremum at 6* = 0, vr and also at 6* = 6*0 for which 



(2) 
vis ■ 



/ ^ {AB-C)x' + Vb^ X + BC 
= {iA-im^2r)x + ^/BC){iA^ 



[m. 



(2)x2^ 



X 



VBC) = 0, 



(117) 



:ii8i 



Then J^^^' 



(119) 



33 



where we have used AB — C = {A — (m[2)^)(^ ~ ("^ii!)^)- Here we will discuss only the 

case with m^2 > ''^12 since the result for m|,2 < ''^12 can be obtained by interchanging 

m^2 m[,2. 

If - C > 0, / > for X > 0, and thus 

— — < for any m^^^ with AB -C >0. (120) 
do 

As jF""*^^' = + m^2 is independent of 9, this means that dJ-'/dO < 0, thus has its 
maximum at 6' = for m|,2, i^x giving AB — C > 0. 

In other case with AB — C < 0, one finds / > for < x ^ Xo? while / < for 
X > Xo, where 



VIS J 



Since 6*0 corresponds to the value of 9 for which x = Xo (for given values of rrL-^}^ and m^), 
this implies 

/ = I - n TV 'f (122) 

[ < for 6^0 < 6^ < TT, ^ ^ 

and thus JF*^*^' decreases as 9 varies from zero to 9o, while it increases as 9 varies from 9o 
to IT. On the other hand, x = a/-B + 4m^ = xo leads to 

,1 _ ^2(2mam, + (mi^;)^ + (ma)^)) 



for which 



-^'^1.=.n=^x + ^12 = -^'^'^'^' (124) 



This means that 9 = 9q corresponds to the boundary between the balanced domain and 
the unbalanced domain. One can show also JF*^^' > jF"'^'^'^^ for 9 > 9q, where the equality 
holds for 9 = 9q, implying that is given by J^^'^^^^ for 9q < 9 < n, while it is given 
by J-"^^^ for < 9 < 9q. Again, combined with that J^^'^^^^ is independent of 9, these 
observations lead to 

— < for AB -C <0, (125) 
o9 

thus has its maximum at 6' = for m^^; ^x giving AB — C < also. 
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B.2 T vs. m|2: 

Let us now examine the dependence of T on vpl^^^^- As we are interested in the events 
giving jF™'^^, we will fix ^ = in the following. We then have 



where 



with 



VIS ) 



The equation (7 = can have the following solutions: 



2A'(A-(mi2)2) + A-(mS)2' 



VIS 



where we have used 

BC - B'C = iA- {m^Ss?) (2A'(A - {m^lf) +A- (m^)') 

and 

'imS^sf - im^sf) {m^Ssf < for ^ m^. 



(126) 



g = (BC - B'C)x^ + VbKJ{B' + l)x + B'BC (127) 



9(m«)2 9(mi2)2 9(m«^2 



If m[2 7^ "^iiL ^6 have also 

B' = {2A' - 1) < 0, 
2A'(A-(ml2)') + A-(mS)' > 0, (130) 

and thus xi is always positive, while the sign of X2 is determined by the sign of A — (m^2)^- 



If m^2 > "^12; ^6 have ^'^^'^i = + m|,2, and thus 

^^unbal 



1 forml2>mi2. (131) 



dm. 



(1) VIS VI 



vts 
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As for the behavior of JF^'^^ since A — > "^12 > ''^12' o'^ly X = Xi can be a 

physical solution of 5^ = 0, for which 



VBC{1 - 2A') 



B + Aml = '—-^ -— . (132) 



VIS J 



This equation is solved by = m^o^, which means g = at = m^o. It is also 
straightforward to find that > for > m^o and g < for < m^o, leading to 

d^^""^ ^ I < for < m^o and > ^^vL (^^33^, 
dm^2 I - for > m^o and > ^vl- 

On the other hand, if < m^o and ^ = 0, is always given by T^^^. We then find from 
f frnl) and (fT33D that 

^ I < for < m^o and m[,2 > ^^^^^ 
ami2 for > m^o and mj,2 > mj,2. 

In other case with m|,2 < "^12' have ^"'^t'^i = + m|,2, so 

^^unbal 

^ (1) = for mj,2 < mSJ. (135) 

In this case, A — (m^2)^ can be either positive or negative. If it is positive, again x = Xi 
is the only solution for g = 0. Then, one can repeat the analysis leading to (11331) . and 
find 

9^ f < for < m^o, mj,2 < ml,2, A - (m|l)2 > ^^^^^ 



dm^ I > for > m^o, < "^12, ^ - ("^12)' > 

However, if A — (m|,2)^ < 0, both x = Xi ^-^id x = X2 become a good solution of 
g = 0, and X2 > Xi- Again, x = xi is obtained if and only if = m^o. As for the other 
solution X = X2, i-e. 



it is straightforward to find that it gives rise to 

^'^1,=,,="^x + "^12=^"'^'^\ (138) 

which means that X = X2 corresponds to the boundary between the balanced domain and 
the unbalanced domain. In Fig. [131 we depict the 2-dimensional event space of {m^^^^, m|,2) 
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for ^ = and > m-j^o, showing the balanced domain (B) and the unbalanced domain 
(UB). 

The above observation implies that JF is always given by J-"^^^ if x < Xi? fo^^ which 
< m^o. One can show that g < in such case, so 

< for < m.o, < m^, A - (m^)' < 0. (139) 



(1) — " ""X ^ ""Xi' 

However, if x ^ Xi? for which > m^o, the function g can have either sign. We find that 
(yf > for Xi < X < while (? < for % > X2- As x = X2 corresponds to the boundary 
between the balanced domain and the unbalanced domain, this means that g > when 
JF is given by T^^^, so 



(1) am^'^ 

(1) / _(2) /I ^^(2)^2 



for > m^o, m^/ < m^^, A - (m^.^) < 0, Xi < X < X2- (140) 
If X ^ X2 with m^2 < ''^iil) is given by jF"'^^^' = ' -^^"^^ 



for > mjo, mi,2 < mj,^, ^ - {rri^^lf < 0, X > X2- (141) 



Combining all of the above results together, and also taking into account that T is 
3 ex( 

dT 



invariant under the exchange of m^2 ^'^'^ "^iiL finally obtain 



dm. 



< for rriy < m^a and any mi*^„ 

- X Xi J VIS (^-^^2) 

> for > m^o and any m^/^. 



=0 

In Fig. (fT^ . we depict the pattern of the 2-d vector field dJ^ / dm~^]^ for both cases with 
< m^o (x < Xi) Ei'iid "^x ^ ^xl ^ ^2)) showing the above result explicitly. 
Together with the observation that has its maximum at 6' = for given values of m'^jg 
and m^, the above result assures that the global maximum of the mT2 function T over 
the 3-dimensional event space parameterized by 6} is given by 

^»"K)=(5: ::";^;:*! (143) 



where 



- J- [m^ig - m„jg , m„j^ - m^^^ , - U, m^j. 



= ^(mi2 = m-r,"^S = <r,^ = 0,mJ. (144) 
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(0 = 0) 




Figure 13: Division of the {fTil'l,m\^^l) into the balanced solution region (B) and 
unbalanced region (UB). 




Figure 14: The gradient of JF for (a) > m^o and 6 = 0, (b) < m^o and 6 = 
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